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II. Number counts and galaxy colours in Bj, R, and K 
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Abstract. We present a medium deep survey carried out in the three filters Bj, R and K. The survey covers 
homogeneously the central square degree around the Northern Ecliptic Pole (NEP) down to a completeness limit 
of 24.25, 23.0 and 17.5 mag in Bj, R and K, respectively. While the near infrared data have been presented in 
the first paper of this series, here we concentrate on the optical data and the results based on the combined 
BjRK -data.. The unique combination of area and depth in our survey allows to perform a variety of investigations 
based on homogeneous material covering more than ten magnitudes in apparent brightness. We analyze the 
number counts for point-like and extended sources in Bj and R to determine the slopes in dlogA^/dm and 
to test for possible breaks therein. While we can confirm the slopes found in previous works with a higher 
statistical significance, the largest uncertainty remaining for the amplitudes is galactic extinction. We determine 
the colour distributions of galaxies in Bj — R and R — K down to Bj — 24.0 and K = 18.0 mag, respectively. 
The distributions in both colours are modeled using galaxy spectral evolution synthesis. We demonstrate that 
the standard models of galaxy evolution are unable to reproduce the steady reddening trend in R — K despite 
flawless fits to the colour distributions in the optical {Bj —R). The BjRK data collected over a large area provides 
the opportunity to select rare objects like candidates for high-redshift galaxies and extremely red objects (EROs, 
R — K > 5.0) and to determine their surface density. Our EROs are selected at an intermediate magnitude 
range and contain contribution from both galactic as well as extragalactic sources. At K < 16.5 mag, where a 
morphological classification is possible, the stellar component dominates the sample. 
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1. Introduction 

Understanding galaxy evolution relies on large and 
homogeneous sets of data. Inhomogeneities introduced 
by stacking surveys with different depth to maintain 
comparable numbers of sources over a large dynamic 
range are an important limitation in testing those models 
of galaxy evolution that are designed to reproduce the 
number counts in different filters. Such models have been 
developed extensively over the last decade (see Koo & 
Kron 1992). Models incorporating luminosity evolution 



Send offprint requests to: M.W. Kummel 

* Based on observations collected at the German-Spanish 
Astronomical Centre, Calar Alto, operated by the Max-Planck- 
Institut fiir Astronomie, Heidelberg, jointly with the Spanish 
National Commission for Astronomy 

** In partial fulfillment of the requirements for a Ph.D., car- 
ried out at Landessternwarte Heidelberg 
*** Email;swagner@lsw.uni-heidelberg.de 
Correspondence to: kuemmel@mpia-hd.mpg.de 



have been found to explain the number counts down 
to the faintest levels (Metcalfe et al. 1995, Pozzetti et 
al. 19961 ). As pointed out by Gardner (1998) however, 
the colour distribution contains more information about 
the state of evolution than the pure number counts, and 
the modeling of the colour distributions is a good test 
for the evolution models which explain number counts 
in individual filters. Up to now such modeling of galaxy 
colours has mostly been done for bright, nearby samples 



(Berlin & Dennefeld 1997) or for deep samples with small 
number statistics (Pozzetti et al. 
al. pOOO| ). 



1996, McCracken et 



We carried out medium-deep surveys in the optical and 
near-infrared regimes. They cover one square degree and 
were performed in the optical Bj- and i?-bands as well 
as with a near-infrared K filter. The 95% completeness 
limits in Bj, R, and K are 24.25 mag, 23.0 mag, and 
17.5 mag, respectively. The near-infrared survey and 
results based exclusively on /iT-data have been presented 
in Kiimmel & Wagner ( |2000| ) (hereafter paper I). The 
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Table 1. The coordinates in the ecHptic, galactic and su- 
pergalactic systems and measures of ISM column density 
towards the NEP 



02000: 18h 00m 00.0s 


S-iooo'- 66° 33' 38.6" 


lii: 96°. 38 


hii: 29°.81 


SGL: 33°. 30 


SGB: 38°. 34 


Hl-column-density^ 


= 3.9-4.7 X lO^Vm'^ 


IR-lOO-^i-cirrus^ 


= 1.8 - 3.0MJy/steradian 




= 0.05 mag 



complete coverage of the iiT-survey by the i?-data, both in 
depth and in the large area, allows us to study the colour 
evolution of our i^-selected as well as optically selected 
sample with a high statistical significance. We report 
on tests of the models proposed for galaxy evolution, 
specifically, we study whether they can reproduce the 
galaxy colours and their variation with brightness. 
One particular field of interest is to determine the 
surface densities of extremely red objects (EROs) in an 
intermediate range of magnitudes. EROs are objects with 
R-K > 5.0 mag (Cimatti et al. |1999i Daddi et al. ^OOOD 
and include galactic and extragalactic populations. While 
the galactic EROs population are late type stars (M6 

|199§) 



or even later, see Leggett 1992 and Wolf et al 



there is no unique explanation for the extragalactic 
EROs. Among the different scenarios (possibly all of 
which contribute on some not yet determined level) are 
galaxies with an old stellar population at high redshift 
with a strong 4000 A break. For a redshift of z > 0.85 this 
break falls between the R- and the K- filter bandpasses, 
resulting in very red colours. Another suggestion for 
EROs are starburst galaxies or active galactic nuclei at a 
redshift 1 < z < 2. In this case reddening by inte rstella r 
dust alters the observed SEDs (Thompson et al. 1999|) . 
A third possibility, which is probably less important 
in the magnitude range covered in our survey are very 
distant quasars, where the Lyman break is redshifted 
to A > 700 nm. In all cases the objects lie at moderate 
to high redshifts and give important clues on galaxy 
evolution and their star formation history. Evidence that 
the largest fraction of the extragalactic component of the 
EROs-population are high-z ellipticals comes from the 
high clustering amplitude suggested in recent surveys 
(Daddi et al. |2000| ). 

The EROs search in our medium deep and medium wide 
survey at the NEP bridges the gap between the large area 



multi colour surveys like DENIS (Epchtein et al. 1997) or 
2MASS (Skrutskie, et al. |199S| ) and the deeper surveys on 
smaller fields, e.g. CADIS (Thompson et al. 1999, Huang 
et al. ^OOlD , or Daddi et al. ( |2000| ). While we can continue 
the DENIS-sea rch fo r low mass stars to larger distances 
(Delfosse et al. 1999 ), we can detect the bright end of the 
extragalactic EROs population, which are identified in 



deep surveys like CADIS. 

The specific field used for our studies is the Northern 
Ecliptic Pole (NEP), which is special in having been 
surveyed intensively by scanning satellites (ROSAT, 
IRAS). Our deep counts shall be used to identify sources 
in deep X-ray and far-infrared surveys (Brinkmann et 
al. 1999, Hacking & Houck 1987) and study their broad 
band energy distribution. 



2. Observations and data analysis 

2.1. Observations 

Two medium deep surveys were carried out with the 3.5-m 
telescope on Calar Alto, Spain in two observing runs from 
July 21-25, 1993 and August 6-8, 1994. During both cam- 
paigns the telescope was equipped with a TEK CCD (CA 
#7) in the prime focus. The CCD has 1024 x 1024 pixels 
with an image scale of 0.403"/pixel. The coordinates of 
the NEP and other field parameters important for survey 
work are given in Table |l|. To cover the central square 
degree around the NEP with the field of view (FOV) of 
6.9 X 6.9arcmin^, an equally spaced grid of 9 x 9 and 
10 X 10 exposures was taken in Bj and R, respectively. 
The exposure times of the individual frames are 10 min 
in Bj and 6 min in R. In order to obtain homogeneous 
photometry over the whole field we carried out a snap- 
shot survey in both bands with short exposure times and 
a large field of view in photometric conditions. For this 
purpose we used the 2.2-m telescope on Calar Alto with 
the focal reducer CAFOS (Meisenheimer 1996). Together 
with the SITe CCD (CA #1) CAFOS has a circular FOV 
with 13arcmin diameter and 0.531"/pixel. In each band 
the snapshot survey was performed on a grid of 4 x 4 ex- 
posures with integration times in Bj and i? of 3 min and 
2 min, respectively. Although the snapshot survey does not 
cover the complete field of the deep exposures, there is 
sufficient overlap with every frame of the deep survey to 
define a common photometric zeropoint. 
All data were obtained using the so-called R6ser-BV 
and R6ser-R2 filters. The R6ser-BV filter {Xcenter = 
497.7 nm and AA — 155.9 nm) is similar to the Bj filter 
199|). The R6ser-R2 (see Roser & 



(see GuUixson et al. 



^ Dickey fc L ockman ( [l990| ) 

^ ISSA ( |199^ ) 

3 Schlegel et al. (|l998|) 



Meisenheimer 1991) avoids the strong OH emission lines 
at A > 760 nm, which contaminate the standard R filters, 
with a sharp cutoff at 740 nm. 

In both surveys the median value of the full width at half 
maximum (FWHM) of the point spread function (PSF) is 
1.5". 

2.2. Standard reduction 

The single raw frames were de-biased and flat-fielded. For 
every observing run a bias frame was constructed using 
frames with an integration time of sec, taken with closed 
CCD-shutter. 

A "super-fiat-field" in each band was obtained by using 
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all exposures taken in that filter. The de-biased exposures 
were normalized and the super-flat-field was computed 
from the median of the data values in every pixel. 
Dark-subtraction could be neglected since none of the 
CCDs displayed significant dark-currents. 



faint sources, respectively (see paper I for a detailed dis- 
cussion). The transition from L total to L;/ca was chosen 
at an object size < 28arcsec^. For those objects the flux 
Lifca was measured in the corresponding aperture of 6" 
diameter. 



2.3. Photometric calibration 

Photometric calibration was obtained by observing several 
standa rd fie lds from Christian et al. ( 1985 ) and Odewahn 
et al. (199S). The _B-magnitudes of the stars listed there 



were transformed to Bj using the equations given by 



Gullixson et al. ( |1995| ). 

We observed standard fields at different airmasses in all 
nights when the snapshot survey was carried out. Instead 
of computing a zeropoint for every exposure of the snap- 
shot survey from the extinction curve individually, we used 
the large overlap between adjacent frames to enhance the 
homogeneity of the photometry. In every overlap region 
bright, unsaturated stars were identified in each pair of 
neighboring exposures to determine the differential zero- 
point between the two exposures. Following the method 
developed by Glazebrook et al. (1994) this system of differ- 
ential zeropoints for every overlap was then transformed to 
a single differential zeropoint for each exposure. In photo- 
metric conditions those differential zeropoints of the expo- 
sures originate from different extinction, hence extinction 
correction is done explicitly. We computed the absolute 
zeropoint of every exposure by adding a constant value, 
which was determined from x^-minimization of the differ- 
ences between the zeropoints computed differentially and 
the zeropoints derived from the extinction curve. 
The zeropoints from the snapshot survey were then trans- 
ferred to each field of the deep survey individually, using 
several stars in each case. 

While we did not find a significant colour term for the 
R6ser-R2 filter, the transformation 



Bj^b + 0.23 * {Bj - R) 



(1) 



was applied to convert the instrumental magnitude b to 



2.4. Object detection and photometry 

Object detection, photometry and morpholo gical c lassifi- 
cation was carried out with FOCAS (Valdes 



1994, Jarvis 



& Tyson 1981). The reliability of the detection process was 
extensively tested on simulated images generated with the 
iraf-package noao.artdata. As in paper I, the threshold pa- 
rameters were chosen such that only < 1% of the objects 
found in the artificial images were not real objects, result- 
ing in a reliability of > 99% for the detected objects. This 
was achieved by setting the FOCAS parameters such that 
after a convolution with the FOCAS built-in digital filter, 
at least nine connected pixels are required to have inten- 
sities > 2.8 cr to be recorded as an object. 
Similar to paper I we used the FOCAS-total flux Lfotai 
and the flux measured in an aperture Li fca for bright and 



2.5. Morphological classification 

The morphological classification into point-like and ex- 
tended sources was done in both the Bj- and the R sur- 
vey. It is based on the FOCAS-resolution classifier (Valdes 
1982 ). This classifier fits a series of templates, which are 



basically derived by scaling the width of the image PSF to 
each object. The scale of the best fitting template is then 
a measure of the resolution of the object and the classifi- 
cation is made from this scale value. 

Like every classifier based on the object shape, the 
FOGAS-resolution classifier does not give reliable results 
for sources near the completeness limit. Towards lower 
signal-to-noise ratios the extended parts of galaxies pro- 
gressively vanish in the background noise. This leads to a 
misclassification of extended sources as point-like objects. 
This misclassification affects sources closer than 1.5 mag 
to the completeness limit. 

Point-like sources are stars, distant quasars, and nuclei 
of galaxies at low and intermediate redshift which have 
steep luminosity profiles such that the width of the nu- 
clear profiles down to the level of sky noise is significantly 
smaller than the PSF. Down to our levels of completeness, 
the surface density of quasars is about 85deg~^, which is 
more than an order of magnitude lower than the density of 
stars accordin g to t he Bahcall-Soneira (Bahc all & Soneira 



1980| , Bahcall |1986| , Mamon & Soneira |1982D model of the 



Galaxy. Nucleated dwarf galaxies, such as M32, would be 
included in our survey out to distances of 200 Mpc. The 
scale length of the nucleus would be 0.2", and the surface 
brightness of the extended emission would be lost in the 
sky noise. Such nucleated dwarf galaxies would be classi- 
fied as point-like sources for 90 % of the volume sampled. 
It is still unlikely, that dwarf galaxies present a significant 
contribution to the point-like sources, since our survey 
only includes one major galaxy out to the distance of 200 
Mpc. The majority of point-like objects that are not stars 
are distant galaxies with small angular size. In the Hubble 
deep fields (HDF, HDFS) 12 % of aU extended objects 
down to i? = 23 mag have scale-length that would render 
them unresolved at our resolution. At brighter magnitudes 
this ratio cannot be determined reliably due to small num- 
ber statistics. 

To derive the number counts of extended objects we used 
a statistical source-classification in the range of unreliable 
FOCAS-classification. We extrapolate the number-counts 
of bright point-like sources by assuming dlogiV/dm = 
const and compute the number of extended objects by 
subtracting the expected number of point-like sources 
from the total number of objects. To test the assumption 
for the statistical classification we computed the number 
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Table 2. The source counts in Bj 



Filter 


mag 






N * 

^ ' star 






14.625 


1.00 


1.00 




1.0 




15.125 


1.00 


1.00 




1.0 




15.625 


1.00 


1.00 




1.0 




16.125 


0.00 


0.00 




1.0 




16.625 


5.02 


2.25 




1.0 




17.125 


8.03 


2.83 




1.0 




17.625 


15.07 


3.89 




1.0 




18.125 


32.15 


5.67 




1.0 




18.625 


34.16 


5.84 


342.58 


1.0 




19.125 


59.27 


7.70 


389.80 


1.0 




19.625 


93.43 


9.66 


453.09 


1.0 




20.125 


161.74 


12.72 


551.54 


1.0 




20.625 


262.21 


16.19 


583.69 


1.0 




21.125 


404.87 


20.12 


693.20 


1.0 




21.625 


727.36 


26.97 


879.05 


1.0 




22.125 


1273.87 


35.70 


970.48 


1.0 




22.625 


2250.38 


47.44 


1075.96 


1.0 




23.125 


3885.19 


62.33 




1.0 




23.625 


7286.75 


85.36 




1.0 




24.000 


10775.98 


146.80 




0.8 




24.250 


15936.32 


178.52 




0.3 



Table 3. The source counts in R 



Filter 


mag 








area 


R 


12.875 


0.97 


0.97 




1.0 




13.375 


0.00 


0.00 




1.0 




13.875 


0.97 


0.97 




1.0 




14.375 


0.00 


0.00 




1.0 




14.875 


1.94 


1.37 




1.0 




15.375 


3.88 


1.94 




1.0 




15.875 


11.65 


3.37 




1.0 




16.375 


22.33 


4.66 


288.70 


1.0 




16.875 


37.88 


6.06 


362.09 


1.0 




17.375 


45.64 


6.66 


393.30 


1.0 




17.875 


62.15 


7.77 


490.41 


1.0 




18.375 


132.07 


13.03 


551.59 


1.0 




18.875 


239.86 


15.26 


624.42 


1.0 




19.375 


390.39 


19.46 


757.47 


1.0 




19.875 


652.59 


25.16 


911.88 


1.0 




20.375 


1028.41 


31.60 


960.43 


1.0 




20.875 


1691.10 


40.52 


1199.05 


1.0 




21.375 


2460.88 


48.88 


1386.00 


1.0 




21.875 


3420.84 


57.63 




1.0 




22.375 


5406.16 


72.45 




1.0 




22.875 


9252.67 


94.78 




0.6 



of stars expected according to the Bahcall-Soneira model 
in Bj and R. Figs. ^ and ^ display the number counts 
of point-hke sources in Bj and R as open circles. Shown 
as a solid hne are the counts expected for the Bahcall- 
Soneira model. There is a good agreement between the 
theoretically expected number of stars and the detected 
number of point-like sources in the range of reliable clas- 
sification. We actually detect more point-like sources than 
expected, confirming a contribution of up to ten percent 
of unresolved galaxies to the list of point-like objects. For 
the stars no significant change in d log TV/dm down to 
Bj = 24.25 mag and R = 23.0 mag is expected in the 
model. This confirms the validity of the assumption in 
the range of unreliable FOCAS-classification and justifies 
our statistical classification. In the following we will treat 
point-sources as stars, despite the above mentioned con- 
tamination. 



2.6. Completeness 

The completeness function / is derived as described in 
paper I. We fit 



Table 4. The source counts in K 



fimag) 



exp 



/ mag — mag^o \ 



1 



(2) 



to the normalized number counts N(m)/N(TO) to deter- 
mine mag^Q, where the number of detected object is half 
the expected one (as inferred from the normalization) . The 
parameter b was found to be independent within the range 
of image quality experienced in our survey at 6 = 0.26, 
while mag^Q varied in both data sets (Bj and R), refiect- 
ing the different image quality for the individual frames. 
We established the relationship between mag^o and the 



Filter 


mag 








area 


K 


7.25 






1.93 


0.9 




7.75 






1.93 


0.9 




8.25 






5.78 


0.9 




8.75 






5.78 


0.9 




9.25 






6.74 


0.9 




9.75 






10.59 


0.9 




10.25 






17.34 


0.9 




10.75 


0.92 


1.02 


23.12 


0.9 




11.25 


0.00 


0.00 


46.23 


0.9 




11.75 


1.93 


1.44 


52.01 


0.9 




12.25 


2.89 


1.76 


86.68 


0.9 




12.75 


2.89 


1.76 


118.47 


0.9 




13.25 


15.41 


4.07 


182.03 


0.9 




13.75 


20.23 


4.66 


275.45 


0.9 




14.25 


34.67 


6.11 


341.91 


0.9 




14.75 


86.68 


9.65 


450.74 


0.9 




15.25 


112.69 


11.01 


549.95 


0.9 




15.75 


265.82 


16.91 


723.31 


0.9 




16.25 


488.31 


22.91 


889.93 


0.9 




16.75 


911.44 


31.65 


1093.54 


0.9 




17.25 


1462.21 


48.96 


1214.90 


0.6 



basic parameters of image quality, i.e. background and 
FWHM of the PSF for the Bj- and the i?-data inde- 
pendently, and used this relationship and the correspond- 
ing zeropoint to calculate magso for every survey image. 
Finally the completeness limit for each image is set at 
magcompi = 'Tiagso — 0.6. At this level the completeness 
function is /{magcompi) > 0.9, and the slope-normalized 
number counts are indistinguishable from 1.0. 
As in paper I the completeness function f(mag) was not 
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Fig. 1. The number counts for extended and point-like 
sources in Bj. The solid line is the expected star counts 
according to the Bahcall-Soneira model 
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Fig. 2. The number counts for extended and point-like 
sources in R. The solid line is the expected star counts 
according to the Bahcall-Soneira model 



used to correct number counts fainter than magcompi- 
However objects down to mag^o were taken to match the 
objects found in different filters and to determine the 
colours of objects. This is justified since those sources 
in the range [magcompi — "^50] which have actually been 
found, form a statistically selected subsample of sources 
detected with the high reliability of 99%, even if not the 
entire population is included. Down to 77150 of the individ- 
ual fields we detect 58 500, 53 900 and 13 000 sources in 
the Bj-, R- and X-survey, respectively. 



the corrections Aa^^ = Aa^; and ASu^^ = J2i ^^ki 
{Aas, = Efc(-Aafci) and ASb, = for B,) 

and applied them to all positions in the whole frame. This 
procedure was iterated once to minimize the contribution 
of a frame with bad astrometry to its overlaps. 
We estimated the accuracy of the astrometry using bright 
stars in the large (~ 1') overlap regions between adjacent 
fields of the i?-survey. The accuracy of the object posi- 
tions was, even for the faintest sources, determined to be 
< 0.8" in each coordinate. 



2.7. Astrometry 

To obtain a precise absolute astrometry for the detected 
objects the Guid e Star Cata log 1 . 2 (GSC) wa s use d 
(Lasker et al. |l990i Russel et al. |l990i Jenkner et al. |l990|) . 
There are 430 GSC-stars in the area covered by our sur- 
veys, and the number of GSC-stars per exposures ranges 
from 1 to 11. All GSC-stars except the brightest one (the 
planetary nebula NGC6543, my = 9.8 mag) were taken 
to establish the astrometry. We identified the GSC-stars 
on the survey images and determined plate constants for 
every observing run and filter using gnomonic projection 
(see Eichhorn 1974). Using the plate constants and the po- 
sition of the GSC-stars we then computed the equatorial 
position of every image center, and, finally, the equatorial 
positions of all objects found on the images. 
Because of the different spacing of the exposures the po- 
sition of an object in Bj and R is based on a different set 
of GSC-stars. The astrometry can therefore considered to 
be independently derived for the Bj- and the i?-survey. 
We took advantage of this in order to compensate for the 
variable number of GSC-stars per image and to enhance 
the overall homogeneity of the astrometry. For every over- 
lap between an individual _Bj -frame B^ and an individual 
i?- frame R/ we computed Aa^; and A5ku the mean offset 
between the Bj- and i?-coordinates of bright stars in right 
ascension and declination, respectively. Then we derived 



2.8. Removal of cosmic ray objects 

Since most of the survey area is covered only by one 
frame, the usual approach to remove cosmics with weight- 
ing maps or masks (Nonino et al. 1999, Arnouts et al. 



1999) could not be followed. Instead we identified sources 
caused by cosmic ray events and removed them from the 
object lists. 

The criterion whether an object is real or just a cosmic 
ray hit is the concentration of the brightness distribution 

cone ~ 2.5 * log((ic/9.0 — ispht) / {ssbr /3.0 j) + mag (3) 

The PSF sets an upper limit in the brightness concentra- 
tion for real objects. Cosmic ray hits usually have higher 
values of eonc, since their shape is not determined by the 
PSF and a cosmic ray event usually affects only a few pix- 
els. Cone is easily derived from the FOCAS-parameters 
Lc (core luminosity), ispht (isophotal brightness) and ssbr 
(sky-noise) together with the object brightness mag (the 
nomenclature of the FOCAS-parameters follows Valdes 
p^82| ). 

Looking at the eone-mag distribution of all objects from 
an individual image the locus of cosmic ray objects could 
be identified easily and the objects could then be removed 
from the lists. 
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Fig. 3. The galaxy counts in Bj from the NEP compared 
to counts from various other surveys 
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Fig. 4. The galaxy counts in R from the NEP compared 
to counts from various other surveys 



3. Object counts 

3.1. Extended and point-like sources in Bj and R 

Figs. |l| and ^ display the object counts in Bj and i?, re- 
spectively. In both figures the counts of point-like sources 
are marked with open circles, the counts for extended 
sources are denoted with the filled symbols. We made no 
attempt to correct the counts beyond the incompletent- 
ness of the individual frames. The counts were only derived 
from fields complete to the specific magnitude. Tables^ 
and 1^ give the number counts (in o6jects/0.5mag/deg ) 
Bj and R, respectively. The counts for extended objects 
are listed in column 3 and counts for point-like objects 
in column 5. The fourth column reflects pure Poissonian 
error of the counts for extended objects. The last column 
gives the area of the sub-survey complete to the specified 
depths. To complete the information concerning number 
counts in our surveys we added the corresponding data for 
the A'-survey in Table ^. 

The point-like sources in Figs. |, I and in Tables |-| are 
only given down to the magnitude of reliable classification. 
The deeper counts of extended objects have been derived 
by subtracting the expected number of point-like objects 
from the counts of all objects (see Sect. 2.5 and below). 
No number densities can be given at the bright end of 
the point-like sources because the objects saturated the 
CCD. The solid line in Figs. |^ and || shows the theoret- 
ically expected stellar counts according to the Bahcall- 
Soneira model (Bahcall & Soneira [ToSOl Bahcall |1986D . 
To calculate the model counts in Bj we transferred the 
original T^-counts with the model B — ^-colours and equa- 
tions given by GuUixson et al. (1995). For the i?-counts we 
changed the code according to Mamon & Soneira ( 19821 ). 
No attempts were made to improve the fit to our data by 
changing the parameters of the model. This is justified by 
the good agreement between our counts and the model. 
The Bahcall-Soneira model does not take into acco unt th e 
so called thick disc introduced by Gilmore & Reid ( 1983|) . 
However, the scope of the comparison done in this paper 



is not to test a particular model of the Galaxy. The agree- 
ment between the counts for point-like objects and the 
Bahcall-Soneira model is taken as confirmation of the sta- 
tistical classification applied to the total c ount s to derive 
the fraction of extended objects (see chap. 2.5). 



3.2. Number densities of galaxies 

In Figs. H and ^ we compare counts of galaxies at the 
NEP with published counts from other surveys in Bj and 
R, respectively. The reference data are from Bertin & 



Dennefeld ( |1997| ), Jon es et a l. (|l99lD , Huang et al. ( |200lD 
and Metcalfe et al. (|1991|, |1995|). AU i?j-counts except 



Huang et al. (2001) were either performed in a i?j-filter or 



transformed to Bj using equations given by the authors. 



To the Huang et al. (2001) photometry we applied the 
transformation Bj — B — 0.19 mag, according to Bertin & 
Dennefeld ( |1997| ) for B-R^ 1.0 mag. In the Figs. | and 
^ the slope 0.5 and 0.4 is subtracted from the logarithm 
of the counts to expand the ordinate and to make differ- 
ences between the counts clearly visible. In order to do a 
quantitative comparison we fitted power-laws of the form 



N{mag) = a * lob^mag-c) 



(4) 



to the data. Tables ^ and ^ show the results of those fits. 
Both Tables give the literature reference and the survey 
area, respectively, (in deg^) in the first two columns. 
The next columns show the slope b, amplitude a and the 
magnitude range in which the fits were done. In Bj and 
R the fits were done for c = 22.0 mag and 20.0 mag, re- 
spectively. The faintest magnitude bins are not included 
in the fit. The numbers in those bins might be affected by 
Eddington-bias (Eddington [1913 ), and could influence the 
fits presented in Tables ^ and |6| significantly since their 
high values are associated with small (relative) errors. 
We fitted different power-laws above and below 19.4 mag 
to the i?-band data from the NEP, since there is a clear 
break in the number counts at this level. While the slope 
is almost 0.5 for the bright magnitudes, it fiattens by more 
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Table 5. A quantitative comparison of galaxy counts in 
different -surveys 



Table 6. A quantitative comparison of galaxy counts in 
different i?-surveys 



survey 


area 


slope 


amplitude 


range 


Bert in 


140 


0.464 ± 0.001 


663 ±5 


16.0-21.0 


Jones 


2.1 


0.442 ± 0.003 


801 ±5 


18.96-23.46 


NEP 


1.0 


0.479 ± 0.005 


566 ±7 


14.38-23.63 


Metc.91 


0.079 


0.491 ± 0.009 


748 ± 27 


18.88-24.38 


Huang 


0.19 


0.473 ± 0.006 


764 ± 18 


16.75-24.75 


Metc.95 


0.005 


0.396 ± 0.001 


1125 ± 143 


22.37-26.87 



than 0.1 towards fainter magnitudes. 
The slopes in Bj are in good agreement with other sur- 
veys with the exception of Metcalfe et al. ( 1995 ). As can 
be seen in Fig. ||, the change in slope at i? ~ 24 (Arnouts 
et al. 1999, Williams et al. 1996) flattens the slope in the 
deep surveys of Metcalfe et al. (1995). While the slope of 



survey 


area 


slope 


amplitude 


range 


Bert in 


140 


0.537 ±0.001 


579 ±3 


14.5-19.5 


NEP, 


1.0 


0.480 ± 0.020 


397 ± 30 


14.13-19.38 


NEPb 


1.0 


0.498 ± 0.044 


424 ± 54 


17.88-19.38 


Jones 


3.0 


0.360 ± 0.002 


448 ±4 


18.13-22.13 


NEPc 


1.0 


0.368 ±0.006 


368 ±5 


18.88-22.38 


Metc.91 


0.079 


0.370 ± 0.008 


432 ± 21 


19.0-23.5 


Huang 


0.19 


0.357 ±0.009 


397 ± 14 


19.25-22.75 


Metc.95 


0.006 


0.399 ±0.017 


315 ±55 


21.75-25.25 



ture of such additional components from number counts 
alone. Surveys in multiple colours provide further insight. 

4.1. Matching of the Bj-, R- and K-catalogs 

The colours of objects are derived by matching the inde- 
pendently produced catalogs in Bj, R and K. For every 
entry in the i?-catalog we searched the i?j-catalog for an 
object closer than 2" to the i?-position. If there exists such 
an entry in Bj, The R and the Bj object are considered 
to be identical, and an entry in the Bj — i?-colour catalog 
is made. The same procedure was applied to the K and R 
catalogs to assemble the i? — i^-colour catalog. Finally, we 
matched the Bj catalog with the colour catalog in R — K 
to get the full colour information of all objects. Since the 
magnitudes derived as outlined in Sect. 2.4 are "total" 



the NEP-counts in Bj is comparable to the slope of other 
surveys to the same limiting magnitudes, the amplitude is 
at a = 566 very low compared to the others, which show 
values around 720. 

In R the NEP-counts clearly resolve the break in the slope 
at 19.0 mag leading from 0.5 at the bright counts to 0.37 
at the faint end. As is the case in Bj the slopes at the NEP 
agree well with the other surveys, but the amplitudes are 
lower. 

The main reason for the low amplitudes in both the Bj- 
and the i?-filters can be attributed to the low galactic 
latitude and therefore high extinction value. While typi- 
cal ex tragalactic survey fields used e.g. in Metcalfe et al. 
( |1995D have Eb-v = 0.02 the extinction at the NEP is 
Eb-v = 0.05. This can be translate d into a fading of 
0.12 mag and 0.07mag (Schmidt-Kaler |l982[ ) of the NEP- 
sources with respect to sources from other surveys in Bj 
and R, respectively. This is supported by the fact that at 
the longest wavelength K there is no such effect (see pa- 
per I) while the difference in the amplitude is strongest in 
Bj, the shortest wavelength. 

Assuming all of the offset is due to extinction, would lead 
us to shift our Bj- and i?-counts by ABj — 0.2 and 
Ai? = 0.15 mag, respectively, corresponding to an extinc- 
tion of /S.Eb-v = 0.05 more than those adopted for the 
fields observed e.g. by Metcalfe et al. (1995). The extinc- 
tion at the NEP then would have to be Eb-v — 0.07, 

0.02 higher than the values from Schlegel et al. (1998) FOCAS-classifier in one filter. For Bj — R we took the 

classification in Bj , and for R — K the classification in R, 



magnitudes, the colours were calculated as the difference 
between the magnitudes in each passband. 
We took 2" as the largest distance for the identification 
of objects in two catalogs. At distances > 2" an increas- 
ing number of object pairs would enter the colour-catalogs 
which are just a positional coincidence hy chance of two 
individual objects in the catalogs of each passband. 



4.2. The colour-magnitude diagram in Bj 

R-K 



R and 



Figs. 1^ and |^ show the colour-magnitude diagrams of ex- 
tended objects in Bj — i? and R — K, respectively. In those 
diagrams the morphological classification is based on the 



(see Tabl e |^), or the extinction towards the Metcalfe et 
al. (|l995| ) fields would have to be neglig ibly small. If the 
extinction given by Schlegel et al. ( 1998 ) is taken into ac- 
count, the difference between the counts at the NEP and 
other surveys reduces to an acceptable amount of ^ 10%. 



4. Galaxy colours 

Changing slopes in number counts can lead to identifica- 
tions of additional components contributing to the entire 
sample of objects, but it is impossible to constrain the na- 



since the PSF of the X-images is grossly undersampled 
(see paper I). Obviously the statistical extension of the 
FOCAS-classification for the number counts (see Sect || 
and paper I) could not be applied to the individual en- 
tries in Figs. ^ and ^ 

Because of the mi sclas sification at the faint end niBj > 
23.5 mag (see Sect. ^.5[) some actually extended objects in 



Bj — R were marked as point-like sources. Those sources 
are not represented in Fig. || and the true population 
is underestimated there. Xn R — K this effect starts at 
R > 22.5 mag and affects only the red faint end in Fig. || 
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Fig. 5. The colour- magnitude diagram of extended 
sources in Bj—R. In the inset median colours and standard 
deviations for magnitude bins indicated by the horizontal 
bars are plotted. 
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Fig. 6. The colour-magnitude diagram of extended 
sources in R—K. In the inset median colours and standard 
deviations for magnitude bins indicated by the horizontal 
bars are plotted. 



to a rather negligible degree. 

As already argued in Sect. 2.£ sources with magnitudes 
down to mag^Q in Bj , R and K enter the colour-magnitude 
diagrams. The sharp cutoffs in Fig. || at the right (faint 
end) and lower right (blue-faint end) reflect the limited 
depth of the Bj- and i?-exposures. Because of the colour- 
term in the transformation of the instrumental- to the 
Bj-magnitudes (see Equ. 0) the cutoff at the faint end is 
not parallel to the ordinate. Due to the large difference in 
the limiting depth of the R- and fT-survey the i?-cutoff in 
Fig ^ is almost unrecognizable in the sparsely populated 
faint red end of the colour-magnitude diagram. 
Objects with unusual colours were individually inspected 
for errors in the reduction or matching process. Whenever 
a problem was detected, e.g. a different splitting of neigh- 
boring objects in the two passbands or a problem in the 
determination of the background near bright stars, the ob- 
jects were excluded from the colour-catalogs. In Bj — R 
and R — K the numbers of removed objects is 30 and 20, 
respectively. Therefore all objects in Figs. || and ||, even 
if isolated in colour, are real objects with good photome- 
try and colours. There are 23 ODD and 4 900 points in the 
colour-magnitude diagrams Fig. || and ^, respectively. 



While the blue objects are assumed to be halo stars sev- 
eral kpc above the galactic disc, the red objects can be 
identified as faint M-dwarfs within the disc, in immedi- 



ate vicinity of the sun (Bahcall 1986, Robin & Crezc 1986, 
Baraffe et al. 199§D . The two populations are not well sep- 
arated but are connected with a less densely populated 
regime at Bj — R ^ 1.7 mag. 

Compared to point-like objects the distribution of ex- 
tended objects in the two-colour diagram is much broader. 
Furthermore, the extended objects usually have a redder 
R — K colour. Following Huang et al. (1997) this allows 
a separation of point-like and extended objects based on 
colours alone in the blue part to Bj — R < 1.6 mag at both 
sides of the line drawn in Fig. |^. In the red part such a 
colour based separation is no longer possible, since the lo- 
cus of the point-like objects then lies above the separating 
line in the colour region populated by extended objects. 
IMost of the 10 % contamination of unresolved galaxies in 
our list of point-like objects are likely to have the same 
distribution in the two-colour diagram as the objects clas- 
sified as extended. They would hence be the dominant 
contribution in the sparsely populated regime with blue 
Bj-R and red R-K colours. 



4.3. The two-colour diagram 

Fig. ^shows the two-colour diagram Bj—R vs. R—K with 
the point-like and extended objects in the left and right 
panels, respectively. As in Figs. ^ and ^ the classification 
into point-like (5 300 objects) and extended sources (4 200 
sources) in Fig. |^ is based on the FOCAS-classification 
in the i?-band. As in the colour-magnitude diagrams, 
only objects brighter than mag^o in Bj, R, and K were 
matched and plotted in Fig. 0. 

Most of the point-like objects in the left part of Fig. are 
concentrated along a well defined line of ^ 0.3 mag width. 



4.4. Colour trends in Bj - R and R-K 

In order to study the colour evolution we derived the me- 
dian colour and its standard deviation in bins of apparent 
magnitude. The values are given in Table |^ and displayed 
in the insets of Figs. || and |^. The horizontal bars mark 
the widths of the bins in magnitude. For the last two bins 
in Fig. 1^ a standard deviation could not be determined, 
since a significant population of blue objects are beyond 
the depth of the i?-survey. In a similar way the last bin in 
Fig. ^ is affected by very red objects being missed in the R- 
survey. Nevertheless it is possible to put those B^-objects 
without i?-counterparts at the blue end of the colour dis- 
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Fig. 7. The two colour diagram for extended (right) and point-like (left) sources The line drawn in both panels divides 
point-like and extended objects in the region Bj — R < 1.9. 

tribution to compute the median displayed in Fig. |^. Table 7. The median colours [medBjR.RK) and standard 

The median Bj — R colour of the galaxies remains con- deviations {o'BjR,RK) in Bj —R and R~ K for magnitudes 



stant at {Bj — R) — 1.36 from the brightest objects down 
to m^. = 22.3 mag. Then follows a rapid evolution to 
bluer colours, reaching {Bj — R) ~ 0.92 at the faintest 
bin. This evolution is triggered by the onset of the pop- 
ulation of faint blue galaxies^ with nriBj > 22.5 mag and 
{Bj ~ R) < Omag. While the full population is present 
down to m^. > 23 mag, only its reddest part can be seen 
in Fig. ^ at fainter magnitudes, because the bluer ones 
have no counterparts in the i?-survey. 
In i? — if there is a steady trend to redder colours towards 
fainter magnitudes. The median of {R — K) = 2.60 mag 
at the bright end changes to {R — K) = 3.82 mag at 
K = 17.8 mag. In the last two bins the evolution to red 
colours seems to level off with the median colour remain- 
ing almost constant. 



slices in Bj and K, respectively 



Bj -range 


medsj R 




Kj -range 


medRK 


(yRK 


16.0-18.0 


1.27 


0.23 


10.5-13.5 


2.60 


0.39 


18.0-19.0 


1.36 


0.35 


13.5-14.5 


2.74 


0.35 


19.0-19.5 


1.31 


0.40 


14.5-15.0 


2.94 


0.60 


19.5-20.0 


1.23 


0.39 


15.0-15.5 


3.22 


0.67 


20.0-20.5 


1.36 


0.45 


15.5-16.0 


3.38 


0.57 


20.5-21.0 


1.34 


0.46 


16.0-16.5 


3.50 


0.66 


21.0-21.5 


1.38 


0.53 


16.5-17.0 


3.66 


0.75 


21.5-22.0 


1.36 


0.59 


17.0-17.5 


3.79 


0.85 


22.0-22.5 


1.36 


0.62 


17.5-18.0 


3.82 


0.87 


22.5-23.0 


1.28 


0.70 








23.0-23.5 


1.11 










23.5-24.0 


0.92 











4.5. An upper limit to Bj dropouts 

An obvious aim for wide-angle surveys is the derivation of 
number densities of rare classes of objects. Cosmologically 
important targets are highly redshifted targets which can 
be identified as drop-out objects when the Lyman edge 



is redshifted to long wavelengths, out of the bandpass of 
individual filters. We determine limits to the surface den- 
sity of candidates for high redshifted objects. This allows 
constraints on the bright end of the luminosity function 
of highly redshifted sources. The if -band limits are not 
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Table 8. The density of high redshift candidates from the 
Lya decrement in the Bj — i?-colour 







area 


cand 


Pcand 


Mas (1350A) 


< 24.375 


< 21.775 


0.3 


131 


451 


-24.00,-25.08 


< 24.125 


< 21.525 


0.8 


312 


371 


-24.25,-25.33 


< 23.875 


< 21.275 


1.0 


287 


287 


-24.50,-25.58 



The densities are given in A''deg 



faint enough to include R — K colour as selection cri- 
terion, and we are confined to the Bj — R index, which 
does not provide a unique identification of highly red- 
shifted sources. Nevertheless, it is interesting for deep 
wide-angle surveys to determine the surface density of 
candidate sources. All i?-band sources above the com- 
pleteness limit of our sample have a counterpart in the 
Bj catalog. The only exceptions to this are a small num- 
ber of faint sources close to very bright stars, where the 
decomposition of faint sources and halo has different ef- 
ficiency in the two bands (see section 4.2). There are no 
true drop-outs in our sample. The amount of the decre- 
ment at the Lyman break is discussed controversially. For 
non-active galaxies without Lyman a lines of high EW, 
a break of about 2. 6 mag 4 ^ magnitudes has been deter- 
mined (Steidel et al. 1999 ). In the filter-system used in our 
survey, the Lyman a-break is between Bj and R for ob- 
jects with z ~ 3.8. Thus all objects with Bj —R > 2.6 mag 
can be regarded as candidates for z > 3.8 galaxies. In 
Table || we give the surface density of those candidates in 
subsurveys with different depth in Bj . The last column of 
Table |^ gives the upper limit in absolute magnitude for 
an object with R = Rum at z = 3.8. The two values re- 
fer to the cosmologies {qo,Xo) = (0.5,0.0) and (0.1,0.0); 
Hq — 50km/sec/Mpc, respectively. Out of the 9464 point- 
like as well as extended sources down to R — 21.0 mag 



we detected no Bj-ha.nd dropout over the entire 1 deg 
field. The reddest objects are at Bj ~ R ^ 3.5 mag. 



According to Steidel et al. ( 1999| ) Lyman-break colours 
G — R > 3.6 mag for redshifts z > 4.15 — 4.45, depend- 
ing on the extinction. Transforming this criterion to our 
filter set, we expect colours Bj — R> 3.6 mag for objects 
with z > 4.35 — 4.65. Using the formulas given in Steidel 



et al.( 1999) we can compute from the limit in appar- 
ent i?-magnitude an upper limit for the absolute magni- 
tude of objects at z ^ 4.5. Depending on cosmology, there 
are no objects with M^ab (1280 A) < -25.0, -26.2 mag for 
(go,Ao) = (0.5,0.0), (0.1, 0.0); i?o = 50km/sec/Mpc re- 
spectively. 



5. Extremely red objects 

5.1. The sample of EROs 

Starting from the colour-magnitude diagram Fig. ^ 
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Fig. 8. The colour-magnitude diagram of the EROs. 
Extended and point-like objects with reliable classifica- 
tion are marked as filled square and triangles, respectively. 
i?-band dropouts are marked with their lower limit. 



those objects, we applied additional selection criteria to 
candidates from Fig. ^. With K < Kcompi we avoid any 
kind of false or spurious detection in the i^-band. Only 
matches with a separation less than 1 arcsecond from K- 
and i?-band objects are accepted for EROs (as compared 
to 2" for Fig. H). This criterion rejects matches between 
one object from close, but resolved, object pairs in R and 
their combined, unresolved counterparts in K. Such con- 
stellations are caused by the large pixel size and resolu- 
tion in K and redden the objects systematically. Since the 
merged -fC-object has a different position with respect to 
both single objects in R, the stronger criterion concern- 
ing separation efficiently removes such mismatches. In R 
we accepted for the EROs every object as counterpart, in 
contrast to Fig. || where only sources with mag < mag^o 
were considered. While this might result in matches with 
non-existing sources in R, no false EROs are produced 
since the R — K colour of a solid detection in K can only 
become bluer. Finally, both authors individually checked 
the ERO candidates for signs of errors in detection, pho- 
tometry and matching of the counterparts in R and K. 
Only objects confirmed by both authors are considered as 
EROs. As i?-dropouts, objects which only have a lower 
limit in their R — i^-colour we considered only objects 
with K < Kcompi- To test for errors the objects were in- 
dividually checked on the if-images. 



5.2. The surface density of EROs 



Fig. g shows the colour-magnitude diagram of the EROs in 
our survey. Unfortunately the threshold for EROs in R—K 
is not very well defi ned, a nd the values vary from R~ K > 
5.0 (Cim atti et al. |l999|) to R - K > 6.0 (Thompson et 
we se- al. 1999 ). Therefore all objects with R — K > 5.0 are in- 



lected the population of extremely red objects or EROs. To 
give a very conservative estimate of the surface density of 



eluded in Fig. ||. ERO objects without detection in R {R- 



dropouts) are included with their lower limit in {R—K) 



lim 



M. Kummel and S. Wagner: Number counts and colours in Bj, R, and K 11 



Table 9. The density of EROs-object as a function of Table 10. Description of the SEDs used to model the 
limits in i^T-magnitude and R ~ K colour. galaxy colours. 



mag 


Pr-k>5 


PR-K>5.5 


PR~K>6.0 


area [deg]^ 


type 


metalicity 


star formation rate 


^form 


evolution 


K < 15.0 


0.09 


0.06 


0.03 


0.93 


El 


2.5 X solar 


exp., T = 1 Gyr 


15 


y 


K < 16.0 


0.21 


0.09 


0.03 


0.93 


E2 


1. Ox solar 


exp., r = 1 Gyr 


15 


y 


K < 17.0 


1.43 (1) 


0.43 (1) 


0.21 (1) 


0.91 


Sa 


1.0 X solar 


exp., r = 4 Gyr 


15 


y 


K < 17.5 


5.46 (7) 


1.91 (4) 


0.91 (3) 


0.61 


Sbc 


0.4 X solar 


exp., r = 7 Gyr 


15 


y 












Scd 


0.2 X solar 


const. 


15 


y 


* The densities are 


jiven in 10 ^arcmin ^. 




Irr 


0.2 X solar 


const. 


1 


n 



(computed via Rcompi ~ K) in Fig. In Table ^ we give 
the surface density of the our EROs for different limits 
in if-magnitude as well as colour R — K. Only EROs in 
sub-surveys complete down to the if-limit specified in the 
first column of Table ^ are taken to compute the surface 
densities. The area of those sub-surveys is given in the 
last column of Table || (see also paper I). The i?-dropouts 
were taken with their respective (i? — K)u,-a to compute 
the surface densities in Table ^. The number of i?-dropouts 
is given in Table ^ within parentheses. 
Only 11 EROs out of the 146 from Fig. ||are bright enough 
to allow a reliable morphological classification in R. While 
the only extended object is marked with a filled square, 
the point-like sources are given as filled triangles. 
The two main contributors to our EROs-population are 
late type stars and galaxies at high redshift [z > 0.8). 
As the available information on morphology suggests, the 
bright e nd is dominated by stellar objects. According to 
Leggett ( 1992 ) a colour R — K > 5.0 mag is expected for 
stellar types M6 and later. With typical absolute mag- 
nitudes of Mk = 9.5 mag and Mk = 11.5 ma g for M6- 
dwarfs and L-dwarfs, respectively (see Leggett 1992 and 
Reid 1999) and we detect these objects out to a distance 
of 400 pc and 160 pc. 

The surface density of extragalactic EROs at the depth of 



our survey is completely unknown. Thompson et al. (1999) 
give a surface density of 0.04arcmin^^ down to K < 
19.0 mag. This is more than 5 times higher than our value 
at R — K > 6.0 for the total population at our limit 
K < 17.5 mag. The high surface density of 0.7arcmin^^ 



as given by Eisenhardt et al. (199£) in their sample down 
to K < 20.1 mag gives clues to a fast decline of the den- 
sity towards brighter magnitudes. Therefore only a few 
out of the 16 objects in the reddest and deepest interval 
of Table ^ might be of extragalactic origin. Deeper studies 
of the red population would require data with both better 
spatial resolution and wavelength coverage. 

6. Modeling of the galaxy colour-distributions in 

Bj - Rand R-K 

The col our trends seen in the colour magnitude diagrams 
in Sect. 4^ are caused by a combination of different stel- 
lar populations sampled in galaxies of higher redshift (and 
therefore observed at fainter magnitudes) and colour cor- 



rections introduced through cosmological effects, the so 
called k-corrections. A method to compute colour 
distributions taking into account those effects in- 
volves the so called pure luminosity evolution models 



(Gardner 1998, Pozzetti et al. 1996). The data in our sur- 
veys give the unique opportunity to compare two observed 
colour distributions with those models. Our large number 
of sources allows us to reject models which do not repro- 
duce the observed distribution with a high significance. 
Note, that 'm R — K the comparison is not affected by an 
incomplete coverage of red if-objects in R. 

6.1. Modeled colour distributions 



This modeling is based on theoretical spectral energy dis- 
tributions (SEDs) computed with e volut ionary synthesis 
techniques by Bruzual & Chariot ( 1993 ). To derive the 
colour-distributions from the input parameters SEDs, 
luminosity function, cosmology and SED-mix (or type 
mix) w e used the program ncmod developed by Gardner 
(1998). The basic characteristics of the SEDs used in terms 
of corresponding galaxy type, metallicity, star formation 
rate, and epoch of first stars Zform are given in Tab. ^ 
The last column indicates whether passive evolution of 
the galaxy luminosity is taken into account or not. In 
all models presented here the cosmological parameters 
Ho = 50km/sec/Mpc, qo = 0.02 and Aq = 0.0 have been 
used. Likewise, all models consider internal absorption by 



dust according to Wang (1991). 
For the theoretical distribution of Bj 
Bj-based luminosity function of Loveday et al. (1992) and 



R we assume the 



a type mix of 0% El, 10% E2, 10% Sa, 15% Sbc, 45% Scd, 
and 20% Irr. 

In R — K we show two mo dels u sing the if-luminosity 
function from Gardner et al. (1997). The type mix is 16% 
El, 16% E2, 28% Sa, 29% Sbc, 5% Scd, and 6% Irr for 
model 1, and 25% El, 25% E2, 36% Sa, 10% Sbc, 3% Scd, 
and 1% Irr for model 2. 



6.2. Comparison of modeled and observed colour 
distribution 

Figs. I and show the comparison between the mod- 
els and the observed colour distribution in Bj — R and 
R — K. The data are shown in solid, the models in 
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Fig. 9. The colour-distribution of galaxies in Bj — R (solid line) in comparison with the distribution of the model 
described in the text (dashed line). 




the dashed and dotted hues. While the observed distri- 
butions have been corrected for galactic absorption by 
A(i? - K) ^ -0.10 mag and A(g, - R) = -0.08 mag 
(see Tab. |l| and Schmidt-Kaler 1982), the theoretical dis- 
tributions in Figs. ^ and 10 were folded with a Gaussian 
of 0.2 mag FWHM in order to mimic photometric errors. 
In the comparison between models and data we had to 
treat the colour distributions in Bj — R and R — K sep- 
arately since none of the models reproduces the colour 
distribution in Bj — R and R — K simultaneously for a 
single luminosity function and type mix. 
Fig. ^ compares the Bj — i?-colour distribution of galaxies 
in bins of different apparent magnitude Bj to the best 
fitting theoretic al c olour distribution (computed as de- 
scribed in Sect 



gether with the two models outlined in Sect. |6.1| . Neither 
model 1 nor model 2 is able to reproduce the observed 
galaxy colours m R — K over the complete range in ap- 
parent magnitude in K . This is caused by the strong evo- 
lution of the whole population to redder colours as already 



demonstrated in Sect. 4.2. While model 1 agrees with the 



observed distribution at the bright end, but is too blue at 
the faint end, model 2 agrees at the faint end but is too 
blue at bright A'- magnitudes. None of the models can fol- 
low the trend of the observed distributions which become 
redder by 1 mag in the range K = 13.5 — 17.5 mag. 



7. Summary 



y). The agreement of the models with ^^^^ performed a medium deep survey in the filters 



the data is very good throughout the magnitude range 
covered. There are small deviations only in the last two 
panels {Bj > 22 mag). In the Bj —22.5 mag bin an agree- 
ment could easily be reached using a broader Gaussian 
filter, which is justified by the larger photometric errors 
at the faint end. The last panel at Bj = 23.5 mag is seri- 
ously affected by incompleteness in both Bj and R, which 
affects the observed colour distribution. 
Fig. nG shows the R — K-co\out distribution of galaxies to- 



Bj, R and K with completeness limits 24.25, 23.0 and 
17.5 mag, respectively in the central square degree around 
the Northern Ecliptic Pole. 

We have derived the object number counts for point-like 
and extended sources in Bj and R and compared our re- 
sults for the extended sources to other galaxy number 
counts found in the literature by fitting a power law to 
the data. In both filters we confirm the slope in dlogN/dm 
found in other, smaller surveys with values around 0.48 
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Fig. 10. Tlie colour-distribution of galaxies in R~ K (solid line) in comparison with the distribution of the models 
described in the text (modell: dashed hne, model2: dotted line). 



and 0.37 in Bj and R, respectively. Differences in the abso- 
lute numbers, represented by the amplitude of the power 
law fits, can largely be attributed to galactic extinction 
and the limited accuracy of reddening corrections. 
We have determined the colour distribution of galaxies in 
Bj — R and R^ K in a large range of 10 magnitudes down 
to Bj = 24.0 and K = 18.0 mag, respectively. In Bj — R 
the median colour remains constant to Bj = 22.3 mag and 
become bluer to fainter levels. This trend to bluer colours 
marks the onset of the so called " faint blue galaxies" (Ellis 
199? ) . In R-K the galaxies become redder from a medium 
colour R - K = 2.6 a,t K = 12 mag to R - K = 5.8 at 
K = 17.8 mag. 

For the filter system used in our survey we have given 
lower limits for the expected colours (in Bj — R) of 
Lyman-break galaxies at z ~ 3.8. We derive the sur- 
face densities for candidates found in our survey. From 
the reddest objects in Bj — R found in our survey we 
have derived the lower limit for Lyman-break galaxies 
at 2 > 4 to Mas (12280 A) > -25.0 mag (for go, Aq, 
ho = 0.5,0.0,0.5). 

We have determined the surface density of red objects 
{R — K > 5.0 mag) down to K = 17.5 mag on the basis of 
our large field of view (0.9 deg^). Since we are unable to de- 



termine the morphology of most of the objects, the surface 
densities include late type stars (M6 and later) as well as 
the extragalactic EROs. At the bright end (K< 16.5 mag), 
which has reliable morphological classification, point-like 
objects dominate the sample. While the surface density of 
objects with R — K>6 declines by a factor of twenty from 



20.0 mag to 19.0 mag (Eisenhardt et al. 1999, Thompson 
et al. 1999, Daddi et al. 2000), we find the surface density 
at 17.5 mag to be reduced only by an additional factor of 
four. This comparison shows that either the steep decline 
in density from X = 20 to 19 mag levels off to if = 17.5 
or that a large fraction of our extremely red objects are 
stars. 

We have compared the colour distribution of galaxies in 
Bj — R and R — K with theoretical colours based on 
spectral evolution synthesis. We have shown that it is 
not possible to find parameters (e.g. luminosity function 
and galaxy type mix) such that both colour distributions 
are reproduced simultaneously. While a type mix domi- 
nated by late type galaxies reproduces the optical colour 
Bj — R, the strong trend to redder colours in R — K with 
increasing magnitudes can not be reproduced by models. 
A model which fits the data aX K = 17 mag is too red at 
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K = 14 mag, a model which fits at the bright end is too 
blue at the faint end. 
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